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SECTION 1.0 


INTRODUCTION 

This report, entitled "Passive Stabilization of the Long Duration Exposure 
Facility" has been prepared for Langley Research Center under contract num- 
ber NAS 1-13440. It presents the results of a seventeen week study on the 
application of the Magnetically Anchored Rate Damper to gravity gradient 
stabilization of the Long Duration Exposure Facility (LDEF) . The study ob- 
jectives were to perform the analyses and simulations required to investigate 
the use of an existing viscous magnetic rate damper for rate stabilizing the 
LDEF spacecraft. The study itself was broken into three main tasks; Linear 
Performance Estimates, Capture and Danger Requirements, and Performance Pre- 
diction. Each of these tasks was performed for two gravity gradient stabil- 
ization configurations; an axisymmetrlc configuration for two-axis (pitch and 
roll) stability; and a non -axisymmetrlc configuration for three-axis stability. 

The report presents the results by stabilization configuration. 



SECTION 2.0 


SUMMARY AND CONCLUSIONS 

2.1 SUMMARY 

The nominal LDEF configuration and the anticipated orbit parameters are 
shown in Table 2.1. Using these parameters, two linear steady state analy- 
ses were performed; one for two axis (pitch and roll) stabilization, and 
one for three axis (pitch, roll and yaw) stabilization. In each of these 
analyses, the effects of orbit eccentricity, solar pressure, aerodynamic 
pressure, magnetic dipole, and the magnetically anchored rate damper were 
evaluated to determine the configuration sensitivity to variations in these 
parameters. The worst case conditions for steady state were Identified, 
and the performance capability calculated. 

Garber Instability bounds (a linear instability associated with gravity grad- 
ient stabilized spacecraft) were evaluated for the range of configurations 
and damping coefficients under consideration. 


The transient damping capabilities of the damper were evaluated for both the 
two and three axis configurations, and the time constant as a function of 
damping coefficient, magnet strength, and spacecraft moment of inertia 
determined. The capture capabilities of the damper were calculated, and the 
results combined with the steady state, transient, and Garber instability 
analyses to select a damper design. 

After completion of the linearized analyses, the selected configurations and 
damper design were siniulated on a large three axis digital computer program. 



Table 2.1 


Nominal Configuration and Orbit Parameters 


Spacecraft weight 

4180 kg 

9200 lbs 

Configuration 

Cylindrical 

Cylindrical 

Diameter 

4.27 m 

14-0 ft. 

Length 

9.15 m 

30.0 ft. 

Moments of Inertia 



Yaw 

19200 kg -m2 

14200.0 slug-ft2 

Roll 

39300 kg-m2 

29000,0 slug-ft2 

Pitch 

39300 kg-m2 

29000.0 slug-ft2 

Products of Inertia 

0 

0.0 

Spacecraft Magnetic 

1000 pole-cm/ 

1000.0 pole-cms/ 

Dipole 

axis 

axis 

Orbit Altitude 

500 km 

270.0 n. miles 

Orbit Inclination 

28.8 degrees 

28.8 degrees 
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The performance of the spacecraft In steady state was simulated for both 
configurations, as was the transient performance. Capture from an initial 
rate of 0.25 deg/sec was also simulated for both configurations. 


2.2 CONCLUSIONS 

The LDEF spacecraft can be gravity gradient to and three axis stabilized 
within the range of parameters considered. 

Table 2.2 shows the linearized performance estimates for the two axis config- 
uration, with associated error breakdown; and the results of the three axis 
digital simulations. The largest sources of attitude error are the magnetic- 
ally anchored rate damper, which is controllable by design; and aerodynamics. 
Solar pressure is small, and as discussed in Section 3.1.2, can never be 
large. The magnetic errors are based upon a 1000 pole-cm magnetic dipole, 
which was the value estimated for the spacecraft. Magnetic errors will, of 
course, increase if the dipole is increased beyond 1000 pole-cm. The effect 
of larger dipoles can be estimated from the errors provided, since magnetic 
errors are linear with magnetic moment up to a value of 30,000 pole-cm (for 
LDEF) . 

The errors caused by aerodynamics depend upon center of pressure/center of 
mass offset, orbit eccentricity, argument of perigee, right ascension, and 
time period with respect to the solar cycle. The errors shown in Table 2,2 
were calculated for a 1978 time period with the worst position of perigee and 
right ascension (section 3.1,5). The effect of center of pressure/center of 
mass offset and orbit eccentricity is considered in the text, and the change 
in performance associated with changes in these parameters can be estimated from 
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Table 2.2 


Two Axis Stabilization Performance 


Linear Error Estimates 


Error Source 



Pitch 

(deg) 

Roll 

(deg) 

Magnetic Rate Damper 

.48 

.65 

Solar Torque 

.02 

.03 

Magnetic Torque 

.06 

.07 

Orbit Eccentricity 

.43 

- 

Aerodynamic Torque 

1.90 

.28 

Total Error (RSS) 

2.03 

.79 


Simulation - Steady State Results * 

Peak Pitch Error 
Peak Roll Error 

Peak Pitch Rate ** 

Peak Roll Rate ** 

Peak Yaw Rate 


1.50 deg 
.88 deg 

.002 deg/sec 
,002 deg/sec 
,001 deg/sec 


Simulation - Transient Results 

Time to Stop Tumbling 150 hours 

Transient Decay Time Constant 43.8 orbits 


* Magnetic dipole not included in simulation because its 
effects are insignificant relative to other disturbances. 


** Relative to orbital 
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Figure 3.15 and Figure 3.25, respectively. 

The capture ability and transient performance values are shown in Table 2.2, 
and the time histories of the spacecraft motion are shown in Figures 5.4 and 

5.5, respectively. For rates below 0.04 deg/sec, the spacecraft will capture 
rightside up (if initially ri^tside up within + 15 deg per axis). For 0.25 
deg/sec, the spacecraft will tumble. As Indicated in Table 2,2 and Figure 

5.5, the two axis configuration did not stop tumbling even after 150 hours. 

The analysis of the performance (Section 5.1.2) indicates that the spacecraft 
will not stop tumbling because its pitch natural frequency (at large ampli- 
tudes) is orbital, and the danger torque, as well as other disturbance tor- 
ques, excite the spacecraft and prevent stabilization. This condition can 
be eliminated only by changing the moments of inertia of the spacecraft to 
change its natural frequency. 

Table 2.3 contains the satoe information as Table 2,2 for the three axis con- 
figuration, which was obtained from the two axis configuration by adding 
182 kg (400 lb) of mass in specific locations. Again, solar torque is small and 
magnetic errors can be scaled linearly with magnetic dipole. The errors 
resulting from an axial offset of the center of pressure/center of mass are 
the same or less than those for the two axis configuration, but the effect 
on yaw of a radial offset is pronounced. A two inch offset has been assumed 
for the error determination. Figure 4.14 can be used to estimate the effect 
on performance of a larger offset. 

The three axis configuration did capture successfully after 76 hours of tumb- 
ling, with the same damper as the one selected for the two axis configuration. 
Its moments of inertia, which were altered by adding mass, prevented the pitch 
resonance condition from occurring. Since the amount of added mass is small, 
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and the difference between the pitch and roll moments of inertia is small 
(less than 826 kg-m^- 609 slug-ft^), a three axis configuration is likely 
to be the nominal configuration because of packaging limitations. 

The damper design selected for both the two and three axis configuration is 
shown in Table 2.4, The design is consistent with existing damper designs. 

The choice of damping coefficient was based upon previous experience, and can 
be increased or decreased if required. 
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Table 2.3 


Three Axis Stabilization Performance 


Linear Error Estimates 


Attitude Errors 


Error Source 

Pitch 

Cdeg) 

Roll 

(deg) 

Yaw 

(deg) 

Magnetic Rate Damper 

.42 

.55 

.89 

Solar Torque 

.01 

.02 

.22 

Magnetic Torque 

.05 

.06 

.43 

Orbit Eccentricity 

.40 

“ 

- 

Aerodynamic Torque 

1.61 

.52 

8.67 

Total Error (RSS) 

1.75 

.68 

9.06 


Simulation - Steady State Results * 

Peak Pitch Error 
Peak Roll Error 
Peak Yaw Error 


1.38 deg 
.62 deg 
9.30 deg 


Peak Pitch Rate ** 
Peak Roll Rate ** 
Peak Yaw Rate 


,002 deg/sec 
.002 deg/ sec 
.002 deg/sec 


Simulation - Transient Results 

Time to Stop Tumbling 
Transient Decay Time Constant 


76 hours 
46.6 orbits 


* Magnetic dipole not included in simulation because its 
effects are insignificant relative to other disturbances. 

** Relative to orbital 
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Table 2.4 


Damper Design 


WELD^^D JOINTS 



Damping Coefficient (Total) 


.692 n.m.sec (.51 Ib-ft-sec) 


Magnet Strength 


225,000 pole-cm 


Suspension 


35.2 dynes/kg (16 dynes/lb) 



SECTION 3.0 


LINEARIZED ANALYSIS - TWO AXIS CONFIGURATION 

The linearized analysis was divided into four parts; Steady State Perfor- 
mance. Garber Instability, Transient and Capture Performance, and Configura- 
tion Optimization. 

3.1 STEADY STATE PERFORMANCE ANALYSIS 

There are five major sources of attitude error for the LDEF gravity gradient 
stabilized spacecraft; the magnetically anchored rate damper, solar torque, 
magnetic torque, orbit eccentricity, and aerodynamic torque. The magnitude 
and nature of the torques, in conjunction with the spacecraft parameters, 
determines the pointing accuracy and capture capability of the spacecraft. 
The following sections describe each of these torques, and their effect 
upon the spacecraft perfortoance. 

3.1.1 MAGNETICALLY ANCHORED RATE DAMPER 

The magnetically anchored rate damper is a GE developed component designed 
to damp large amplitude oscillations of gravity gradient stabilized space- 
craft. The viscous fluid version of the damper is shown in Figure 3.1. It 
consists of an inner sphere, which contains a permanent magnet and bellows; 
and an outer sphere, constructed of pyrolitic graphite (for diamagnetic 
centering) and aluminum. The space between the spheres, and inside the 
bellows, is filled with Silicone oil, with a viscosity selected to provide 
the required damping coefficient. 

The mechanism of damping depends upon the relative rate of rotation of the 
inner sphere of the damper, and the outer sphere; which is rigidly attached 
to the spacecraft. During spacecraft acquisition and capture, the relative 
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motion is largely the result of motions of the spacecraft. After the space- 
craft has stabilized, however, relative motion continues to exist because the 
magnet follows the Earth's magnetic field, not the local vertical. The damp- 
ing torques then become disturbance torques to the spacecraft and contribute 
to the overall pointing inaccuracy. The amplitude of this error is a linear 
function of the damping coefficient, and decreasing the damping level will 
decrease the contribution to steady state error. However, decreasing the 
damping coefficient will also increase the time to decay to the steady-state 
condition from initial conditions. As a consequence, a tradeoff exists bet- 
ween steady state accuracy and dancing time constant. Section 3.4 discusses 
this tradeoff. 

The steady state performance for the LDEF was determined both from a linear 
dynamic analysis, and a non-linear torque analysis. It is known from these 
analyses and many simulations, that the damper produces sinusoidal torques 
with a wide frequency spectrum, including high harmonics of orbital frequency. 
Consequently dynamic oscillations at several frequencies appear on the space- 
craft attitude from the damper alone.- To determine accurately the steady 
state characteristics, the torques at the frequencies of highest importance 
must be calculated. 

The steady state disturbance torques induced by the magnetically anchored 
rate damper were calculated by an analysis of the earth's magnetic field, 
assuming the spacecraft was perfectly oriented. For sinqalicity, the earth's 
magnetic field was assumed to be a single dipole field oriented along the 
spin axis of the earth, with the damper magnet following the field exactly. 

The rate of change of the magnetic field vector with respect to the orbiting 
coordinate system was determined, and instantaneous torque calculated. The 
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torques were then calculated at equal Increments around the orbit, and fre- 
quency components (i.e,, orbital frequency and the harmonics) computed by 
a Fourier Series analysis. Three torque coefficients were calculated for 
each axis; static, orbital, and twice orbital. Higher harmonics were not 
calculated since they produce negligible attitude errors. 

For a two axis stabilization performance analysis, only the pitch and roll 
coefficients are of importance. The yaw axis has no preferred orientation 
and simply "drifts" under the influence of the damper. The attitude errors 
resulting from the damper pitch and roll torques depend upon the damping 
coefficient, the orbit inclination, and the sensitivity of the roll and yaw 
axes to the disturbance frequencies. The pitch axis is disturbed most with 
torque components at all frequencies. Roll, however, has only orbital and 
twice orbital frequency coefficients, and yaw has only orbital frequency 
torques. The axes and Euler angles are shown in Figure 3.2. 

The relationship between the steady-state damper Induced error and damping 
coefficient Is shown in Figure 3.3 for the LDEF two axis configuration. The 
steady-state errors, especially the pitch errors, are quite low. The worst 
roll error is slightly over one degree for the highest damping coefficient 
considered. Similarly the highest pitch error is 0.5 degrees. This is well 
below the 3-deg pitch bias error at which Garber instability can occur for 
the range of damping coefficients considered, (see Section 3.2). Note that 
the damper does not develop constant yaw torque and hence does not produce 
"infinite" yaw error (or continuous motion). 

3,1.2 SOLAR TORQUE 

Solar torque is the result of the solar force vector, caused by the pressure 
of the sunlight (approximately 47.8 x lO'^ n/n? = iqf/(,lb/ft2) , not passing 

j. 

through the spacecraft center of mass. The location ofi^the’^pej^tfir.-of (mass 
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is well defined, but the effective point of application of the solar force 
(i.e., the center of pressure) Is not invariant, but moves as a function of 
sun angle. For a spacecraft as large as the LDEF, solar torque will be a 
strong function of sun angle, since the solar force vector is close to the 
surface of the spacecraft and the center of mass is on or near the axis of 
symmetry. 

To generalize the study, solar torque was calculated for center of mass/ 
center of geometry differences between 0 and ,762 m (2.5 feet) along the 
axis of symmetry. This range more than encompasses the range of parameters 
likely to be encountered. For exavsxple, if the top half of the spacecraft were 
purely reflective (specular), and the lower half were purely absorptive, 
the effective center of pressure/center of mass misalignment would be 
.366 m (1,2 feet), less than half the maximum value considered. 

The results of the solar torque analysis are shown in Figure 3.4. For this 
evaluation, the sun was assumed to be in the orbit plane, which places the 
torque primarily on the pitch axis. This is the worst sun orientation since 
"out of plane" torque will affect the roll axis, which is stronger than the 
pitch axi^ and would produce less local vertical pointing error. It is 
evident from Figure 3.4 that the largest solar torque error is an orbital 
frequency error on pitch, and that it is negligible (.24 degrees worst case 
considered) . Solar torque does not produce any significant pitch constant 
torque. 

3,1.3 MAGNETIC DIPOLE TORQUE 

In addition to indirectly causing a spacecraft torque through the magnetically 
anchored rate damper, the earth's magnetic field can produce a direct magnetic 


- 16 - 



ftZSi 




••• 


i • ' 

•*1 



:u.:. :.j 


IT.' 

'!:j 








r . ' 

::!.l 

i;:ii 

"ill 








‘i;: 









m 

Or'-i 



« » • • 1 • • » 

ill iill 





i;; .jji: 






w 

ii 

' 1 

. . 1 4..1 . 


i-iiiuii 

•/C-, 

't!*Hr*r- 


SRW 

/r>- 




L'-'-U'’' 






torque. Any magnetic dipole, caused by magnetic materials mounted on the 
spacecraft, will attempt to align itself to the earth’s magnetic field and 
will torque the spacecraft. The ma^itude of the torque is proportional 
to the strength of the magnetic dipole moment, the orientation of the dipole 
within the spacecraft, and the location of the spacecraft with respect to 
the earth. The magnetic field strength decreases with the cube of the orbit 
radius (exactly the same as gravity gradient) , and is twice as strong at 
the poles (north and south) as at the equator. Hence, a spacecraft in a 
high inclination orbit has more magnetic torque on it than a spacecraft 
in a low inclination orbit, and similarly a spacecraft at a high altitude 
has less torque on it than a spacecraft at low altitude. 

To quantitatively determine the magnitude of the attitude errors, the fre- 
quency of the magnetic disturbance torques must be known. At low altitudes, 
the local magnetic field changes as the spacecraft moves in orbit, and the 
magnetic torque changes as a function of time. For orbit periods short 
compared to earth’s rotational period, the torques (pitch, roll and yaw) 
are largely constant (zero frequency) and/or sinusoidal at orbital frequency. 


A magnetic analysis was performed with the same magnetic field assuir5>tions 
as those of Section 3.1.1 except direct magnetic torques were considered. 
Torques appear on all axes, being exclusively sinusoidal at orbital frequency 
on pitch, and constant plus sinusoidal at orbital frequency on roll and yaw. 
In low inclinations the roll and yaw magnetic dipoles are the most effective, 
producing static torques on yaw and roll (respectively), with the pitch axis 
dipole producing only small sinusoidal roll and yaw torques. For a two axis 
gravity gradient spacecraft, a constant yaw magnetic torque implies a pre- 
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ferred yaw orientation, since the dipole will align Itself to the magnetic 
field about the yaw axis. The torcjues on pitch reduce with orbit inclina- 
tion, but in fact never go to zero since the earth’s magnetic equator is not 
coincident with the earth's geographic equator. 

The magnitude of the magnetic torques are directly proportional to the size 
of the spacecraft magnetic moment. In general, this dipole is not known in 
advance of the spacecraft construction, and varies with payload. The magnetic 
dipole of the LDEF spacecraft is not expected to be large because of the nature 
of its payload. Hence, only a single dipole magnitude, 1000 pole-cm per 
axis, was evaluated. The results of the evaluation are shown In Figure 3.4. 
Since magnetic torques are pure torque, not force, the error is independent 
of center of mass position. 

3.1,4 ORBIT ECCENTRICITY TORQUES 

One of the characteristics of a circular orbit is the constant rate of rota- 
tion of the radius vector (a vector from the center of the earth to the orbit- 
ing body). For a spherical earth the radius vector is parallel to the 
local vertical, and a gravity gradient spacecraft will align itself with the 
local vertical and acquire the average rate of rotation. 

The rate of rotation of the radius vector Is not uniform for an eccentric 
orbit, however, but varies from a minimum at apogee to a maximum at perigee. 

The spacecraft will acquire the average rate of rotation of the eccentric 
orbit, but cannot respond to the variations in rate. Consequently errors will 
develop between the axis of minimum moment of inertia and the local vertical. 
The spacecraft will therefore be torqued sinusoidally by gravity gradient, with 
peak torque at apogee and perigee. Since the disturbance torque is a gravity 
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gradient torque, absolute moments of inertia are irrelevant, and the attitude 
error is dependent only upon moment of inertia relationships, and orbit 
eccentricity 

The attitude error resulting from orbit eccentricity is only a pitch error. 
Figure 3,5 shows the attitude error as a function of orbit eccentricity for 
the LDEF. Neither the orbit altitude, nor the orbit inclination, affect 
eccentricity errors. 

3.1.5 AERODYNAMICS 

At altitudes of 500 km (270 nm) the earth's atmosphere is a major source of 
disturbance to the LDEF. The dynamic pressure (the familiar 1/2 pv^) caused 
by the spacecraft's passage through the ratified atmosphere can exceed solar 
pressure by a factor of ten or more. Aerodynamic pressure is directly pro- 
portional to the aerodynamic density, and to estimate the density, a model 
of the earth's atmosphere is required. There are several atmospheric models 
in existence, the most widely used being that of JACCHIA^^\ This model 
defines the atmospheric density as a function of altitude, solar sun spot 
cycle (F 10.7 cm solar flux) and diurnal bulge. The diurnal bulge is a 
"thickening" of the earth's atmosphere due to the earth temperature increase 
associated with solar heating. Typically, the greatest density occurs at 
2 p.m, local time, where the density may be a factor of three higher than 
on the opposite side of the earth (at this altitude). The density differ- 
ence is greater at higher altitudes. > 

Aerod 3 mamlc torques, like solar torques, are dependent upon the spacecraft 
configuration, but are particularly sensitive to orbit altitude, sun spot 
cycle, orbit eccentricity, argument of perigee, and position of the orbit 
with respect to the diurnal bulge. 
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A general aerodynamic study was performed for LDEF, since it was anticipated 
that aerodynamic effects would be the largest source of disturbance torque. 

Four factors considered in the analysis were spacecraft configuration, orbit 
eccentricity, arguii^nt of perigee, and ascending node. The F10.7 centimeter 
solar flux index was selected to be 150, which is approximately 95 percentile 
for the 1978 time period^ and the nominal altitude was 500 km (270 n. miles). 
A drag coefficient of 2.0 was assumed. 

3. 1.5.1 Spacecraft Configuration 

The LDEF spacecraft is symmetrical, and the aerodynamic torques would be 
negligible if the center of aerodynamic pressure coincided with the center 
of mass. The effect of angle of attack between the aerodynamic stream and 
the spacecraft, is insignificant because the angle of attack is relatively 
constant. Any shift in the center of mass from the center of geometry will, 
however, produce an aerodynamic torque, and hence an attitude error. 

The effects on pitch and roll errors of the aerodynamic torques caused by an 
axial offset of the center of mass/center of pressure are shown in Figures 
3.6 through 3.13. Axial offsets were considered the most important for the 
two axis configuration. The torques were calculated for a circular orbit 
with an ascending node position of zero degrees (at autumnal equinox) . 

Figures 3.6, 3.7 and 3.8 show the constant, orbital frequency, and twice orbit- 
al frequency pitch errors resulting from aerodynamic torques on the pitch 
axis. The largest errors are those due to the orbital frequency component 
of the torque, reaching a peak of approximately 5 degrees at a center of 
mass offset of .762 m (2.5 feet). The constant torque error reaches a peak 


- 22 - 






1 



,0 ; AO 2.0 3-C 


I . ’ 

I--;:; 

LUiiXu-;.:. 


^ z>/sr'A/</c^^ or c.m. rroM cr/vr^r, ft. 

. . -. 




!':•■ 



JiB 




i!i:i 

ii.: 

l-VT 

B;:l 







,; 
■ ■ « i 

i. : 



; • 


w / 

. : : 

’~C & 

i:!;' 





• • I • ' 

■ ;:1:t 


rL ; 


; ■ : 



ii-' 



•rZT*. : 

; ■ : : 





’. .-n 

• • 

ni; 

;:;• 





:;;i 


tfif* 

hir 








ili- 









r.'--.) 



J 



" i;:-’ 




Si 







J+tr 

r 1 ; 




. , . . 1 . . . <• 

ii'iilin: 


?:i!!;!i: iS 

Sjfe'il!?;: 

. : ; ; i ; • ■ 





j;:;- 

■:r. 

11 



rUi:.; 

■;!■!■ 


frr 


i-;; 

m 

SI ' 




-S: 















ttBU- 



ajm 




s»j 




















DlSrAMCr* DF C;M.L^^Ftefijfsf^nCE'N'rER’ 


F iGOKE: 3. & 


EFFEcr, Of SHlf T :QF 
PtrCM EiRROa -te :T0. 

TFF. ri. i 

.fiTd'H 

. 1- : i-'-, ; 
















of approximately 2.5 degrees at the same center of mass offset. The value 
of the constant torque error is particularly important with regard to Gar- 
ber Instability (Section 3,2). The pitch error due to the twice orbital 
component of the aerodynamic torque is small (.65 deg. at .762 m) . 


The motions of the roll and yaw axes are coupled, even in the linear model, 
and hence roll errors are obtained from both roll and yaw disturbance torques. 
Figure 3.9, 3.10 and 3.11 show the effects of the axial shift of the center 
mass on the constant, orbital frequency and the twice orbital frequency com- 
ponents, respectively, of the roll error caused by the roll torque. The 
primary cause of the roll torque is the rotational rate of the earth, which 
causes the atmosphere at an altitude of 500 km to have a velocity in inertial 
space; a component of which is along the pitch axis. The resulting torque Is 
therefore along the roll axis. All of the roll errors are small. 

Figure 3.12 and 3,13 demonstrate the Influence on the orbital and the twice 
orbital frequency components of the roll error caused by the yaw torques. 

The constant components are zero because the roll yaw coupling is zero. All 
these errors are negligible. 

For completeness, center of mass/center of pressure shifts along the cylinder 
radius (normal to the axis) were considered. For the two axis configuration, 
such a shift would cause the spacecraft to weathervane about yaw, and assume 
a preferred orientation. Pitch and roll torques would be reduced in this 
manner, but some pitch and roll errors will still result. Figure 3.14 shows 
the results of the analysis, and indicates no significant pitch or roll error 
increase (an axial offset of 0.152 m - six inches - is included in the figure). 
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Figure 3.15 shows the total aerodynamic error in each axis as 


a function of 


axial offset. 

3,1. 5.2 Orbit Eccentricity 

Figure 3.16, 3.17 and 3.18 show the effect of the orbit eccentricity on the 
constant, orbital, and twice orbital frequency components of the aerodynamic 
pitch error. The center of pressure/center of mass offset was selected as 
.152 m (six inches) for this evaluation. This would require a difference in 
drag coefficient of approximately 13 percent between the top half of the 
cylinder and the bottom half of the cylinder, or a center of mass shift of 
.152 m (six inches), or a combination of both. For a cylinder, a 13 percent 
change in drag coefficient is virtually lnq>ossible, hence a center of mass 
shift Is assumed. It is evident from the figures that the effect of orbit 
eccentricity is severe at large (0.02) eccentricities. It must be recog- 
nized, however, that at an eccentricity of 0.0729, perigee is at the surface 
of the earth, well below any reasonable perigee. 

The roll errors caused by aerodynamics in an eccentric orbit are shown In 
Figures 3.19, 3.20 and 3.21; and the roll errors due to yaw torques at 
orbital and twice orbital frequency are shown In Figure 3.22 and Figure 3,23. 
These errors are negligiblei 

The effect of a radial center of mass shift was evaluated for an orbit with 
a 0.01 eccentricity. Figure 3.24 shows the pitch and roll errors assuming 
a .152 m (six inch) axial offset. As in Figure 3.14, there is no pitch error 
increase, but the roll error is much larger than in Figure 3.19 through Figure 
3.23. Figure 3.25 shows the total effect of orbit eccentricity on each axis. 

As Indicated on the previous Figures, other error sources are negligibly affect- 
ed by orbit eccentricity. 
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3. 1.5. 3 Argument of Perigee / Ascending Node 

Figures 3.26, 3.27 and 3.28 show the effect on the pitch error of orbit 
position relative to the sun. Because of its large error, only pitch has 
been considered. These error estiioates were calculated at autumnal equinox, 
which places the diurnal bulge on the equator. Therefore, by changing the 
ascending node position, the effect of both ascending node and argument of 
perigee (best case and worst case positions) are determined. 

Three eccentricities were considered: 0, 0,01 and 0.02, The center of mass/ 
center of pressure offset was again selected to be .152 m (six inches). 

The differences between the worst case and best case performance is low 
at zero eccentricity (approximately .21 deg.), but increases rapidly with 
increasing eccentricity reaching a difference of 4.7 deg. at an eccentricity 
of 0.02, For values of eccentricity anticipated for the shuttle (Section 
3.4), the effects of ascending node and argument of perigee on spacecraft 
performance are small. 

3.2 GARBER INSTABILITY 

All gravity gradient spacecraft have a first order instability which was 
identified by T. Garber of Rand.^*^^ This instability is caused by a con- 
stant torque on the pitch axis, which produces a constant pitch error. This 
error causes the roll axis to go unstable, typically with a long period 
exponential. The instability appears in the linearized equations of motion 
when the linearization is performed about bias pitch, roll and yaw positions. 
Only pitch biases produce the instability, however. 

The pitch bias position at which the spacecraft goes unstable is a function 
of configuration and damping. Increasing the damping will increase the pitch 
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error at which the instability occurs. Only two disturbances cause constant 
pitch errors; the magnetically anchored rate damper and aerodynamics. The 
damper causes a pitch bias, but also supplies damping, and within the limits 
of the current spacecraft designs, produces no instability. The major source 
of constant torque is therefore aerodynamics (Figures 3.6 and 3.16). 

Table 3.1 shows the pitch errors at which Garber instability occurs as a 
function of damping coefficient. At the lowest level of damping coefficient, 
a three degree pitch bias is required. The required pitch bias increases 
with damping coefficient up to a damping level of 1.08 n.m.sec (0.8 Ib-ft- 
sec) , at which point it drops abruptly. The reason for this drop, as well 
as a discussion as to its validity, is given in appendix A, Based upon the 
estimate of the aerodynamic torque, however, Garber instability will not 
occur. 


3.3 TRANSIENT AND CAPTURE PERFORMANCE 
3.3.1 TRANSIENT PERFORMANCE 

To determine the damping and steady state characteristics of the damper, 
the equations of motion of the damper and the spacecraft must be solved simul- 
taneously, Linearization of the equations is the normal analytical approach, 
but the magnetic field orientation changes as a function of orbit position 
and linearization is only approximate (Appendix A), Consequently, a dynamic 
pitch planar analysis is used for most damping calculations, with several 
simplifying assumptions about the magnetic field model. With this approach, 
a simple relationship between damping constant, magnet strength, spacecraft 
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Tatle 3.1 


Effects of Damping on the 
Garber Stability Boundaries 

Magnet Strength = 225000 pole-cms 


Damping Coefficient 

Pitch Bias 

n-m- sec 

Ib-ft-sec 

deg 

0,542 

0.4 

3 

0,678 

0.5 

4 

0.813 

0.6 

5 

0.950 

0.7 

6 

1.08 

0,8 

4 * 

1.22 

0.9 

7 


* See Appendix A 
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moment of inertia ratio, and damping time constant can be determined* The 


relationship requires only two non-dimensional terms: b and M > 

Izz^o 

where b is the damping coefficient, Izz is spacecraft moment of 
inertia, ^o is orbital rate, and M Is the damper magnetic moment (in 
pole-cm). Selection of these terms defines the damping time constant in 
terms of orbits, which remains the same irrespective of altitude. 

Figure 3.29 shows the relationship between pitch damping time and daiiper 
characteristics for the LDEF. The curve represents the pitch performance 
of a spacecraft in a polar orbit about a simple magnetic dipole. More 
Intricate linear analyses, as well as digital simulations, have indicated 
the analysis is a good approximation for high inclinations, and gives a 
rough estimate of the performance at the lower inclinations. 

Three magnetic strengths for the damper are shown in this figure, A magnet 
strength of 225,000 pole-cm is compatible with an existing GE damper design. 
The figure indicates that for damping coefficients higher than 2.71 n.m.sec 
(2.0 Ib-ft-sec/rad) , this magnet is pulled off the Earth's field, and the 
damping performance degrades drastically. Larger magnets improve the damp- 
ing, but only at high damping coefficients, which produce higher damping. 


These figures were developed assuming the damper has no moment of inertia. 
Analyses were also performed assuming that the damper had the anticipated 
moments of inertia. The analyses, contained in appendix A, indicate similar 
performance, except for instabilities on the roll damping curve. 
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3.3.2 CAPTURE PERFORMANCE 

It is evident that as long as the magnet stays anchored to the Earth's mg- 
netic field, it will damp out any oscillations that the spacecraft is subject- 
ed to (see Section 5,1.2, however). But the damper Toagnet may be pulled away 
from the magnetic field by initial angular rates which exceed its design 
capability. The maximum angular deviations of the magnet from the Earth's 
field are inversely proportional to the magnet strength and directly propor- 
tional to the damping coefficient and the spacecraft rates. Past experience 
has demonstrated that the magnet strength should be such that the damper mag- 
net does not deviate by more than approximately 40.0 degrees from the Earth's 
magnetic field (a conservative value). 

Figure 3.30 shows the maximum allowable spacecraft angular rates as a func- 
tion of the damping coefficient for a given magnet strength and maximum 
angular deviation of the magnet. It indicates that for a magnet strength 
of 225,000 pole-cms and a spacecraft angular rate of 0,25 degs/sec, a damp- 
ing coefficient of 0.882 n.m, sec (0.65 Ib-ft-sec) is required to produce 
an angular deviation of 40 degrees. 

An additional design constraint for the magnet arises from the magnitudes of 
the decentering forces it has to tolerate. The existing GE designs can 
provide 35,2 dynes of centering force per kilogram (16 dynes/pound) of the 
magnet mass. Figure 3.31 shows the centrifugal decentering force to be 
withstood by the damper magnet on the LDEF spacecraft. It shows that for 
maximum. angular rates of 0.25 degs/sec the decentering force is less than 
3,52 dynes/kg (1.6 dynes/lb), even if the damper is offset from the center 
of mass of LDEF by as much as 1.97 m (six feet). Thus, the centrifugal 
forces are well under the limit of 35.2 dynes/kg (16 dynes/lb). 
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3.4 CONFIGURATION OPTIMIZATION 

Optimization of the two axis LDEF configuration consists of identifying as 
closely as possible the magnitudes and frequencies of the external torques, 
estimating the spacecraft performance, and selecting the damper characteristics 
required to damp the spacecraft. The damper selected must also prevent the 
Garber instability from occurring, and induce only small spacecraft oscilla- 
tions. The disturbance torques and their effect on the spacecraft perfor- 
mance are contained in Section 3.1, on a parametric basis. It is necessary 
to define, within the range of parameters considered, the conditions the 
spacecraft will encounter. Aerodynamics is the largest source of error for 
the LDEF; depending upon the actual value of center of pressure/center of 
mass offset and orbit eccentricity. An offset of 0.152 m (6 inches) was 
selected as the nominal value, as discussed in Section 3.1.5. The eccen- 
tricity of the orbit depends upon the value that can be achieved by the shuttle. 
The shuttle specification^^\ indicates an insertion capability of + 4 nm 
on the orbit, which at 500 Km (270 nm) produces an eccentricity of .00082, 

Taking a conservative value of 0.002, the eccentricity error is .43 deg., 
and the aerodynamic pitch error (worst case) is 1.9 degrees. The solar 
torque error, with the same offset, is less than ,01 degrees in pitch, and 
less than 0.1 degrees in roll. The magnetic error (for a 1000 pole-cm mag- 
netic dipole in all axes) is ,06 degrees In pitch, and .07 degrees in roll. 


The damping coefficient for the magnetically anchored rate damper was selected 
based upon the desired damping time, Garber instability and steady state 
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error. Consideration was also given to the required damper characteristics 
for the three axis configuration (Section 4.0), It is evident from Figure 
3.30, that a dipole of 225,000 pole-cm Is adequate for the anticipated worst 
case acquisition conditions. Since this value is compatible with existing 
damper designs, the magnet strength was accepted. 

The daii^jing coefficient was selected based upon a tradeoff between damping 
time constant and steady state attitude error, factoring in the Garber in- 
stability restrictions. The limitations on the selection of damping coeffic- 
ient provided by the Garber instability are shown in Table 3.1. Since the 
constant pitch error caused by aerodynamics is 0.5 degrees, and the lowest 
instability angle is three degrees, instability will not occur for any of 
the damping coefficients considered. Hence, only the damping time/steady 
state tradeoff is important. Previous experience with low altitude space- 
craft (NRL, GEOS, etc.) has indicated a damping time constant on the order 
of twenty orbits provides acceptable damping performance with reasonable 
damper induced errors. Consequently, a damping coefficient of 0,692 n.m.sec 
(0.51 Ib-ft-sec) was selected, providing a linear damping time constant of 
22.5 orbits and steady state pitch and roll errors of 0.39 degrees and 0,63 
degrees, respectively. The error budget for the configuration which com- 
bines all of these errors, is shown in Table 3.2. The errors are arrived 
at by root-sum-squaring all the error sources with the same frequency com- 
ponents, and then adding "across" the frequencies. The approach is conser- 
vative and produces somewhat larger errors than would actually be encountered. 
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Table 3.2 


Linear Error Breakdown 
Two Axis Configuration 
Spacecraft Parameters 

Orbit Eccentricity = 0.002 

Position of the Ascending Node and Perigee = 2:00 P.M. 

C.M, Offset Along Pitch Axis = 1.0 ins. 

C.M, Offset Along Roll Axis = 1.0 ins. 

C.M, Offset Along Yaw Axis = 6.0 ins. 

Damping Coefficient = 0.51 Ib-ft-sec/rad 

Spacecraft Magnetics = 1000.0 pole-cms/axis 


Torque 

Error 

Axis 

1 

Orbit Freq. 

Error Magnitude Deg, 

Frequency 

Components 

Total 

Orbit 






Eccen- 

Pitch 


1 

.432 

.432 

trlcity 









0 

.521 



Pitch 


1 

1.270 

1.902 




2 

.111 


Aerody- 






namic 









0 

.008 



Roll 


1 

.258 

.283 




2 

,017 

- ■ - - - - - --- - j 
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Table 3.2 (contM) 


Torque 

Error 

Axis 

Error Freq, 

Orbit Freq, 

Error Magnitude Deg, 

Frequency 

Components 

Total 



0 

0.113 



Pitch 

1 


0.274 

.483 



2 


-0.096 


Damper 








0 




Roll 

1 



.648 



2 






0 


0.0 



Pitch 

1 


0.005 

.019 



2 


-0.014 


Solar 








0 





Roll 

1 


0,0X3 

.026 



2 






0 


0.0 



Pitch 

1 


0.060 

.06 



2 


0.0 


Magnet 








0 


0.021 



Roll 

1 


0.048 

.069 



2 


0.0 



Total Pitch Error 2.03 degrees 

Total Roll Error .79 degrees 
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SECTION 4.0 


LINEARIZED ANALYSIS - THREE AXIS CONFIGURATION 

The linearized analysis for the three axis configuration was divided into the 
same four parts as the two axis configuration analyses; Steady State Perfor- 
mance, Garber Instability, Transient and Capture Performance, and Config- 
uration Optimization, 

4.1 STEADY STATE PERFORMANCE ANALYSES 

To provide three axis attitude control, the spacecraft must have three diff- 
erent moments of inertia; the largest on the pitch axis, the smallest on 
the yaw axis, and the intermediate on the roll axis. Attitude stabiliza- 
tion is provided directly by gravity gradient on the pitch and roll axes, and 
by "dynamic torques" on the yaw axis. The "torques" result from the fact 
that the spacecraft prefers to rotate about its maximum moment of Inertia, 
and will move toward that orientation when forced to rotate at orbital rate 
by gravity gradient torques. 

To obtain a three axis control configuration for the LDEF, mass was added to 
the spacecraft in the roll-yaw plane. Four equal masses were placed sym- 
metrically at the points farthest from the pitch axis (Figure 4.1), Since 
the three axis configuration has the same envelope and shape as the two 
axis configuration, the magnitude and frequency of the disturbance torques 
are virtually unchanged. The change in performance from the two axis config- 
uration is due almost entirely to added mass, resulting in increased moment 
of inertia (Figure 4.2) and altered frequency response. The largest effect 
is on the yaw response to constant torques, which is reduced from the infinite 
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response of the two axis configuration. The yaw response to orbital fre- 
quency is only slightly improved, since the large yaw moment of inertia 
of the LDEF already provides good frequency response. 

As indicated in the previous paragraph, the disturbance torques affecting 
the three axis configuration are the same as those affecting the two axis 
configuration. It is therefore only necessary to evaluate the errors pro- 
duced on the three axis configuration by these torques. The largest source 
of torque is aerodynamic pressure, as indicated in Section 3.1.5, and the 
largest contributer to aerodynamic torque is a center of mass/center of 
pressure offset. The effect of an axial offset of the center of mass is 
the same or less for the three axis configuration as the two axis config- 
uration. Figures 4.3, 4,4 and 4.5 show the pitch, roll and yaw performance 
of the spacecraft with an axial offset of .152 m (six inches) ; the value 
selected as nominal in Section 3.4. The effect of a radial center of pres- 
sure /center of mass offset is, however, much more severe on a three axis 
configuration than a two axis configuration, primarily because of yaw. 

The errors are largest on yaw, which is the weakest axis. Its restoring 
torque, for 182 kg (400 lb) of added mass is 174 dyne -cm/ degree compared 
to 15,200 dyne-cm/degree for pitch. 

The effect of a center of pressure/center of mass offset along a radius 
(normal to the centerline axis) is shown in Figures 4.6 through 4.9 as a 
function of added mass. For convenience, all the errors for each axis 
were added, as opposed to root sum squaring. The most significant improve- 
ment is the reduction in "constant yaw error" with increasing mass. The 
primary purpose of adding the mass is to provide a yaw restoring torque. 
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This restoring torque is directly proportional to the difference between the 
pitch and roll moments of inertia. Since the geometry of the spacecraft is 
fixed, and the masses are placed in the optimum position, the difference 
between pitch and roll is directly proportional to the mass. Hence, the yaw 
gravity gradient torque is linear with added weight, and the yaw error is 
inverse with added weight. This effect is evident In Figures 4.6 through 
4.9, where the yaw error decreases directly as the added weight increases. 
With 182 kg (400 pounds) of added mass, for example, a .152 m (six inch) 
center of pressure/center of mass offset produces a constant yaw error of 
approximately 10.1 degrees. At 908 kg (2000 pounds), the error is only 
1.94 degrees. When the added weight contributes significantly to the LDEF 
moment of inertia, a decrease in the yaw sinusoidal errors appears (pitch 
and roll, also) . 


4.2 GARBER INSTABILITY 

Each of the configurations analyzed in Section 4.1 was evaluated for Garber 
instabilities. The results are shown in Table 4.1. It is evident from this 
table that the Increased moments of inertia reduce the effectiveness of the 
damping, and lowers the pitch bias at which the instability occurs. The 
result was anticipated since increasing the spacecraft moment of inertia 
without increasing the damping coefficient, reduces the effective damping, 
and consequently lowers the instability point. The lowest allowable pitch 
bias is three degrees for any level of damping considered. The relationship 
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between the Garber stability boundaries and the spacecraft performance is 
summarized in Figure 4.10 which indicates that Garber instability is not 
a problem. 

4.3 TRANSIENT AND CAPTURE PERFORMANCE 

4.3.1 transient performance 

The transient analysis performed for the two axis configuration was also 
performed for the three axis configuration. The results are shown in Figure 
4.11, 4.12 and 4.13. The general relationship between dan 5 >ing coeffic 

tent and damping time constant established in Section 3.3 is retained, but 
the curves are shifted, showing reduced damping with increased mass. Here 
again, the reduced damping is the result of increased moment of inertia 
for the same damping coefficient. The effect of increased magnet size is 
the same for the three axis configuration as the two axis configuration. 

4.3.2 CAPTURE PERFORMANCE 

The spacecraft configuration is irrelevant to the damper as far as capture 
is concerned (see Section 5.1.2, however). The spacecraft will ultimately 
capture if the magnet is anchored to the magnetic field. Hence, the analy- 
sis performed in Section 3.3.2 for the two axis configuration is directly 
applicable to the three axis configuration. 

4.4 CONFIGURATION OPTIMIZATION 

Optimization of the three axis LDEF configuration consists of selecting the 
amount of mass which must be added to provide yaw stabilization; and select- 
ing the damper characteristics to provide appropriate damping and stability. 
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An evaluation of the disturbance torques acting on the two axis configuration 
(Section 3.4) led to the selection of a 0.152 m (six inch) center of pressure/ 
center of mass offset, for a 0.002 orbit eccentricity value. These values 
have only a small inpact in yaw, however, as indicated in Figure 4,5. The 
principal source of yaw error is a center of pressure/center of mass offset 
normal to the cylinder axis, along a radius. The worst position is along 
the spacecraft pitch axis (Figure 3.2), which maximizes the yaw torque due 
to aerodynamic pressure. For offsets of .101 m (four inches) or less, 
reasonable performance can be obtained by adding 182 kg (400 lb) to 364 kg 
(800 lb) of additional mass. From discussions with Langley Research Center, 
an offset of .025 m (one inch) can be achieved by weight balancing. To 
obtain a 0.025 m (one inch) offset by coefficient of drag differences, how- 
ever, would require a temperature difference of several hundred degrees 
between the left and right side of the spacecraft, or a completely differ- 
ent shape on the two sides. Neither of these conditions is likely, but to 
be conservative, an additional inch was added to the center of mass location 
to produce a two inch center of pressure/center of mass offset. With this 
offset, 182 kg (400 pounds) of additional mass will provide yaw control to 
within nine degrees. The variation of aerodynamic error with center of 
pressure/center of mass offset along a radius is shown in Figure 4.14. 

The selection of damping coefficient was relatively straightforward in view 
of the small change in moment of inertia between the two axis configuration 
and the three axis configuration (with 182 kg). With the same damping coef- 
ficient as the two axis configuration (0.692 n.m.sec), the linear damping 
time constant is 24 orbits. The damper characteristics for the two axis config- 
uration was therefore selected for the three axis configuration. The error 
breakdown for the three axis configuration is shown in Table 4.2. 
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Table 4.2 


Linear Error Budgets 
Three Axis Configuration 

Spacecraft Parameters 


Orbit Eccentricity = 0.002 

Position of the Ascending Node and Perigee = 2:00 P.M, 

C.M. Offset Along Pitch Axis = 2.0 ins. 

C.M. Offset Along Yaw Axis = 6.0 ins. 

Damping Coefficient = 0.51 Ib-f t-sec/rad 

Spacecraft Magnetics = 1000.0 pole-cms/axis 


Torque 

Error 

Axis 

Error Freq, 

Orbit Freq 

Error Magnitude, Deg, 

Frequency 

Components 

Total 

Orbit 





Eccen- 

Pitch 

1 

.40 

.40 

tricity 







0 




Pitch 

1 


1.61 



2 

.10 


Aerody- 


0 

.01 


namic 

Roll 

1 

.43 

.52 



2 

.08 




0 

7.64 



Yaw 

1 

.87 

8.67 



2 

.16 
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Table 4.2 (contM) 



Torque 

Error 

Axis 

cr 

<u 

u 

w 

0 

u 

u 

w 

♦ 

O' 

u 

p* 

•H 

u 

o 

Error Magnitude, Deg. 

Frequency 

Components 

Total 



0 

.10 

.42 


Pitch 

1 

.23 

.42 



2 

.09 




0 

.0 



Roll 

1 

.53 

.55 

Danq>er 


2 

.02 




0 

.0 



Yaw 

1 

.89 

00 



2 

.0 




0 

.0 



Pitch 

1 

.0 

.01 



2 

.01 




0 

.0 



Roll 

1 

.01 

.02 

Solar 


2 

.01 




0 

.19 



Yaw 

1 


.22 



2 

,01 




0 

.0 



Pitch 

1 

.05 

.05 



2 

.0 




0 

.0 



Roll 

1 

.04 

.06 



2 

.02 


Magnet 







0 

.34 



Yaw 

1 

.07 

.43 



2 

.02 




































Table 4.2 (cont’d) 


Total Pitch Error 
Total Roll Error 
Total Yaw Error 


1.75 deg 
.68 deg 
9.06 deg 
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A comment about the moment of inertia distribution for the three axis config- 
uration appears warranted at this point. The basic LDEF spacecraft has been 
assumed to have exactly 39,300 kg-m^ (29,000 slug-ft^) in the pitch and roll 
axes. With the addition of 182 kg (400 lb), the difference between the 
pitch and roll moments of inertia reaches only 825 kg-m^ (609 slug-ft^). 

This value is small compared to the basic spacecraft, and it is estimated 
that the value falls below the packaging accuracy capability. Consequently, 
the spacecraft is likely to have three different moments of inertia, the 
differences of which exceed the value necessary to provide good yaw control. 
If the moments of Inertia can be controlled to be favorable (pitch a max- 
imum, yaw a minimum), and the moment of difference between pitch and roll 
equal to greater than the value calculated for 182 kg of added mass, the 
added mass is unnecessary. The exact moment of inertia difference is not 
critical to performance, and if the pitch moment of inertia is approximately 
the same as the two axis configuration, damping will not be significantly 
affected. 
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SECTION 5.0 


PERFORMANCE PREDICTION 

Performance Prediction consists of simulating the LDEF configuration using 
a three axis digital computer program. The purpose of the simulation is 
to verify the conclusions drawn from the linear analyses. Capture is also 
simulated, because of its non-linear dynamics. 

The computer program used to make the simulations is a large digital computer 
program which simulates the dynamics of the spacecraft in three axes using 
Euler djmamical equations. The attitude of the spacecraft is specified by 
Euler parameters. Solar torque, aerodynamic torque, magnetic torque and 
damper torque are included in the simulation for both circular and eccen- 
tric orbits. The torques are calculated as a function of spacecraft atti- 
tude and orbit position. 

5.1 TWO AXIS CONFIGURATION 

5.1.1 STEADY STATE PERFORMANCE 

Several computer simulations were made to determine the behavior of the 
spacecraft once steady state is achieved. The first simulation included 
only the magnetically anchored rate dair^ser and orbit eccentricity, and was 
selected to illustrate the behavior of the magnetically anchored rate damper 
in steady state. Figure 5.1 shows the last fifty eight hours of the simu- 
lation. Only the end of the simulation is shown to eliminate transients 
caused by errors in initial conditions. The damping coefficient in the 
simulation is .697 n.ra.sec (.51 Ib-ft-sec) . The orbit eccentricity is 0.002, 
with perigee located at the base of the diurnal bulge. 
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The error breakdown. Table 3.2, indicates the pitch axis will show a pitch 
bias of 0,1 degrees and an oscillation of 0,61 degrees. Figure 5.1 indi- 
cates a bias of 0,25 degrees, with an oscillation of 0.63 degrees. There 
is a long period of oscillation on the spacecraft pitch axis, with a period 
of approxiioately twelve hours. This is the effect of the rotation of the 
earth on the magnetic field. 

The behavior of the spacecraft in roll is very different from the linear 
estimate. The largest difference is a pronounced bias of approximately 0.5 
degrees, not indicated by the linear analysts. The roll bias and the yaw 
rotation (obvious in Figure 5.1) are directly connected. If the roll bias 
were the cause of the yaw rate, a roll bias of 0.5 degrees would produce a 
yaw rate of 5.55 x 10'^ deg/sec (orbital rate times the sine of 0.5 degrees); 
approximately 180 degrees in 90 hours. The yaw plot indicates a rate of 
1,23 X 10'3 deg/ sec (180 degrees in 40,5 hours). Considering that the yaw 
rate is causing a roll bias, a yaw rate of 1.23 x 10"^ deg/sec will cause 
a coupling torque of 6.66 x 10“^ lb- ft, which will induce a roll error of 
0.534 degrees. It is apparent that It Is the rate about yaw which is caus- 
ing the roll bias, rather than the reverse. Figure 5.2 confirms the analy- 
sis since the bias disappeared when yaw stopped rotating. The cause of the 
yaw rotation appears to be the result of a linear instability caused by the 
damper (appendix A). The instability is in yaw, and manifests itself as 
a steady yaw rate. The effect on steady state is not serious, since the 
overall effect of the instability is to cause the spacecraft to assume a 
biased attitude (and rate). The roll error (exclusive of the bias) is 
0.63 degrees, in good agreement with the linear estimate. 


C 
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Figure 5.2 shows the steady state perfomance of the two axis configuration 
with all error sources (solar torque, aerodynamic torque, damper, and orbit 
eccentricity) except magnetic torques. The error breakdown for this case 
indicates a pitch error of 2.03 degrees and a roll error of 0.79 degrees. 

The roll error shown in Figure 5.2 is approximately 0.8 degrees (note the 
bias has disappeared), in good agreement with the linearized estimate. The 
pitch axis, however, has a peak error of approximately 1.5 degrees, nearly 
0.5 degrees below the anticipated value. This reduced error is primarily 
the result of reduced oscillation amplitude. The bias has increased approx- 
imately 0.4 degrees between Figures 5.1 and 5.2, compared to the linear 
estimate of 0.52 degrees. 

A significant difference between Figures 5.1 and 5.2 is the position assumed 
by the yaw axis. The spacecraft is weathervaning as a result of a center of 
mass offset of one inch along each of the roll and pitch axes. The weather- 
vaning eliminated the observed effect of the yaw instability, but it cannot 
be established that the instability is eliminated. The instability may, for 
example, cause a yaw bias from the weathervaning null. The overall agree- 
ment with the linearized estimate, however, is reasonably good, since, as 
indicated earlier, the linearized estimate is conservative. Additional 
steady state runs were unnecessary because of the excellent performance, 
even in the worst case, and the close correlation between the linearized 
analysis and the simulation. However, an additional simulation was made using 
a damping coefficient of 1.38 n.m.sec (1,0 Ib-ft-sec) to determine the effect 
of a higher damping coefficient. Figure 5.3 shows the results, and indicates 
a peak error of 1.9 degrees in the pitch axis, and 1.2 degrees in the roll 
axis. This is increase of .4 degrees in pitch, and ,3 degrees in roll due 
to the damper, in agreement with the linearized analysis (Figure 3,3), 
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5.1.2 transient and CAPTURE PERFORMANCE 

The transient and capture simulations were performed with two sets of initial 
conditions. The first set placed 0.04 deg/sec rate and fifteen degrees 
attitude error on all axes. The conditions were selected (from Langley 
Inputs) to verify that the spacecraft would stabilize rightside up. The 
conditions also provided an excellent run for the calculation of damping 
time constant. 

The results of the simulation are shown in Figure 5.4. The spacecraft does 
capture upright, reaching a peak pointing error of 59 degrees at five hours. 

The local vertical at the end of the simulation is 11.3 degrees, indicating 
a transient decay time of 90 hours. This Is not in agreement with the linear- 
ized estimate (23 orbits), and Is directly attributable to the behavior of 
pitch. The time constant for pitch is 43 orbits, nearly twice the anticipated 
value. In general, pitch damping falls off with decreasing orbit inclina- 
tion, and some increase in damping time is normally anticipated. Correlation 
with other simulations had indicated that to an inclination of 45 degrees, 
the reduction is minor. However, the majority of the previous simulations 
were for two axis configurations with pitch and roll moments of inertia much 
larger than yaw (factor of 100). A brief comparison with the results obtain- 
ed from a study performed for Marshal Spaceflight Center indicated closer 
correlation with LDEF simulation results. Other simulations performed during 
this study (Section 5.2.2) indicated similar results. Hence, the linear time 
constant is 43.8 orbits. 

The second set of initial condltons was a capture condition with initial rates 
of 0.25 deg/ sec on all axes, and 15 degrees of attitude error. The rates 
are the maximum separation rates anticipated for the Remote Manipulator for 
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the space shuttle. Figure 5.5 shows the results of the simulation, and 
Indicates that the spacecraft tumbled for 150 hours. Increasing the damp- 
ing from .692 n.m.sec (0.51 Ib-ft-sec) to 1.36 n.m.sec (1.0 Ib-ft-sec) did 
not improve the damping (Figure 5.5). It Is evident from Figure 5.6, that 
the spacecraft Is not damping after the first 60 hours. An analysis of the 
results Indicated that the primary cause of the prolonged tumble Is a space- 
craft resonant condition associated with large amplitude oscillations. Since 
gravity gradient torques are a function of the angle between the local ver- 
tical and the minimum moment of Inertia, a gravity gradient spacecraft will 
act as a spring mass system. At low amplitude oscillations, the natural 
frequency of the spacecraft In pitch for the two axis configuration Is 1.23 
times orbital. At large amplitudes, however, the gravity gradient "spring 
constant" reduces, and the natural frequency of the spacecraft falls below 
the low amplitude natural frequency. At oscillation amplitudes of ninety 
degrees, the configuration's natural frequency Is nearly orbital (the fre- 
quency can be estimated from Figure 5.5). As a consequence, the config- 
uration Is extremely sensitive to torques which disturb the spacecraft at 
orbital frequency, and the damper (as well as aerodynamics) drives the space- 
craft and prevents It from settling. The problem can be eliminated by alter- 
ing the natural frequency of the spacecraft through appropriate moment of 
Inertia change. Altering the damping coefficient, however. Is not effect- 
ive since the level of damping is too low to significantly Impact the 
natural frequency. 

5.2 THREE AXIS COWFIGURATIQH 

5.2.1 STEADY STATE PERFORMANCE 

The three axis configuration was simulated in a manner similar to the two 
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axis configuration. The first simulation included only the magnetically 
anchored rate damper (damping coefficient of .692 n.m.sec (.51 Ib-ft-sec) 
and orbit eccentricity (.002). Figure 5.7 shows the last fifty-eight hours 
of the simulation, and indicates pitch, roll and yaw errors of 0.82 deg. 

0.47 deg, and 3.3 degrees, respectively. The error breakdown indicates 
pitch, roll and yaw errors of 0.65 deg, .55 deg, and .89 deg for this case. 

The pitch and roll estimates are in reasonable agreement, but the yaw axis 
simulation results are significantly different from the linearized estimate. 
The difference in the yaw performance Is primarily the result of a low 
frequency oscillation (twelve hours) associated with the rotation of the 
earth. The once orbital frequency, superimposed on the long cycle, has 
an amplitude of 0.65 degrees, which is in good agreement with the linearized 

estimate . 

The performance of the three axis configuration with all disturbances, is 
shown in Figure 5.8. The peak pitch, roll and yaw errors are 0.94 deg. 

0.47 deg, and 8.2 deg. The linearized error budget Indicates 1,75 deg, 

0.68 deg, and 9.06 deg for the same conditions. The better agreement in 
yaw between the linearized estimate and the simulation is probably the result 
of aerodynamic torque suppressing the low frequency yaw disturbance. The 
low pitch error appears to be Che result of phasing of the disturbances 
since there is an obvious difference in behavior between Figures 5.7 and 
5.8, but little increase in oscillation amplitude. 

As with the two axis configuration, a steady state run was made with a damp- 
ing coefficient of 1.36 n.m.sec (1,0 lb-£t-sec). The results, shown in 
Figure 5.9, indicate pitch, roll and yaw errors of 1.7 deg, 1.2 deg, and 
10.5 degrees. The long period oscillation on the yaw axis is more pro- 
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nounced than in Figure 5.7 and Figure 5.8. 


5.2.2 TRANSIENT AND CAPTURE PERFORMANCE 

The transient behavior of the spacecraft from initial rates of 0.04 deg/ sec 
on all axes and fifteen degrees of attitude error is shown in Figure 5.10. 

The spacecraft captured rlghtside up, and damped with a time constant of 
46.6 orbits; in good agreement with the two axis simulation (Figure 5.5), 

The natural period in yaw is approximately eight hours as shown on the Figure. 
This oscillation is not to be confused with the 12 hour period shown in 
Figure 5.7. 

A second transient simulation was made without external disturbance torques 
to determine the effect of aerodynamics, solar torque, etc. on the 
damping. The spacecraft behavior shown in Figure 5.11, is different than 
that shown in Figure 5.10, but the time constant is basically the same, 

A third transient simulation was made with a damping coefficient of 1.36 
n.m.sec (1.0 Ib-ft-sec), and the results are shown in Figure 5.12. As with 
the lighter damping, the spacecraft captured rightside up; but the time 
constant is 27.7 orbits. 

A capture run was made with 0.25 deg/sec on all axes, the results of which 
are shown in Figure 5.13. Tlie spacecraft tumbled in pitch for approximately 
76 hours, after which it stabilized rightside up and began damping. The pro- 
longed tumble observed in the two axis configuration is not in evidence. The 
natural frequency of the three axis configuration is 0.93 orbital at large 
an5)litudes, and hence the configuration is not in resonance with orbital 
frequency disturbances. As indicated in the Figure, the spacecraft stabilized 
in the rightside-up position, but rightside up capture is problematical with 
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APPENDIX - A 


DAMPER MAGNET DYNAMICS 



NOMENCLATURE 


La3 , 
L&3 

b 

LC3 

LBl 


State matrices; Eqs. (33), (34). 

Damping coefficient matrix; Eqs. (5), (6). 
Scalar damping coefficient; Eq. (6). 

State matrix; Eq. (35) 

Coordinate transformation matrix; Eqs. (1) , 

( 8 ). 


e. 




'} j 


j ” 1 , 2 ) 3 - 


Elements of QcJ » Eq. (2). 


^ Principal moments of inertia of the 

^ ^ spacecraft in the x;;j2 coordinates; 

Eq. (3). 





Principal moments of inertia of the 
magnet in the coordinates. 




KOLL 




/o/tcM 


,r. 




Roll, pitch and yaw moments of Inertia 
of the spacecraft in the three-axis 
configuration . 


tzj 


Identity matrix. 


Universal gravitational constant multi- 
plied by the mass of the Earth. 


M 


Strength of the magnet in pole-cms. 


L03 


Null matrix. 




Orbital frame of coordinates; Eq. (1) , 


•Co 


Geocentric distance of the spacecraft. 


A-1 




; 1= f -4. 

Coefficient matrices of the magnet 
equations of motion; Eqs=> (28), (29). 


; J, »c - l.i.3. 

Elements of ; Eq. (29). 

T 


Maximum yaw stiffness of a magnet; 
Eqs. (25), (30). 

Zis 


Damper torque on the spacecraft; Eq. (5). 



Gravity gradient torque vector on the 
spacecraft; Eqs, (4), (7). 


T T” 

Magnet stiffnesses; Eq. (23). 



Magnetic torque vector on the magnet; 
Eqs. (12), (15), (25), 


, l-,-4. 

Coefficient matrices of the spacecraft 
equations of motion; Eqs, (26), (27). 


) rjflk. — 

Elements of 

f 


Time from the ascending node. 

a 


State vector; Eq. (31), 



Body fixed coordinate system of the 
spacecraft; Eq. (1). 

o( 


Local latitude. 



Orbi t inc 1 Ina t ion . 

# 

1. 


Orbital rate. 

^p , ^ 

r, 

Pitch, roll and yaw angles of the spacecraft 


A-2 



^f>s i ^rs * 

Bias values of • 


Perturbed values of ^ 

from the bias values. 




Spacecraft rate vector in y ^ Z - 
coordinates . 


Damper rate vector referred to the 
axes fixed on the spacecraft. 


Spacecraft rate vector referred to 
the axes fixed on the damper. 

SUPERSCRIPTS 


r-y 

Denotes variables which relate to the 
damper magnet. 

SUBSCRIPTS 


X j, y ^ 

Denotes the elements of a vector in the 
corresponding axes. 

OPERATORS 


c 

Transpose 

• 

( ; 

Derivative with respect to time. 
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EULER’S EQUATIONS OF MOTION 

A coordinate system Is aasiimed to be attached to the orbiting space- 

craft. Let Y^<jf be the orbital coordinate system rotating at the orbital 
rate of the spacecraft about the ^^-axis. The relationship between the 
coordinate systems is defined as 
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where 

the elements 
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= 

S>-i^ ^ Sy, 


^f3 

= 



^2/ 

= 

Cea Sly. &Y ^ 



— 

C.&0 0^ C.^ ^ ^ ^ 




Ce-3 Oy. 



=r 

sx^&y 



s 

S Sy (2&^ ^ ^ 'S 


^23 

rr. 



( 2 ) 
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The angles and are the pitch, the roll and the yaw angular 

displacements, respectively, of the spacecraft and define the rotation of 
the x^2 coordinates with respect to the coordinates. This 

angular relationship between the coordinate systems is shown in Figure Al. 

’pjie ^ -axis is taken along the local vertical and the ^ -axis is parallel 
to the trajectory in the direction of motion of the center of mass of the 
spacecraft. The Z -axis is normal to the orbit plane forming a right 
handed coordinate frame. With this convention, and assuming that the 
system is coincident with the axes of principal moment of inertia of the 
spacecraft, these moments of inertia are given by 

~ 3 ^3 ” ‘ ( 3 ) 

With the preceding assun^tions, the gravity gradient torque vector, , 

in the coordinate system on the spacecraft is given by 

^ ^ ^3 J 

7}, = ^ Ul. --Tj e,, 

where the subscripts ^ 2 refer to the cooresponding axes, and 

= universal gravitational constant multiplied by the mass 
of the Earth, 

a geocentric distance of the satellite. 
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Let cOj be the angular velocity vector of the danger magnet expressed 

in the xy2 coordinate system fixed on the spacecraft. Then the damping 

torque vector T, on the spacecraft created by its motion relative to 

di& 

that of the damper magnet is given by 

4 . = 


where 

« the angular velocity vector of the spacecraft in the xyZ 
coordinates* and 



the damping coefficient matrix for the damper assembly. 


Due to the isotropy and symmetry of the damper assembly, the matrix 
is a scalar matrix, and is of the form 


= bii:i 


( 6 ) 


where b scalar damping coefficient and 

Qx] is the identity matrix. 

Torques on the spacecraft from all other sources are neglected as the effects 
of only the damper are under investigation. Hence, the equations of motion 
of the spacecraft are given by 


I. 


-h 


) CO^ 


b ( 







) 

+ 

w (Wj 

" ^c\sy.) 

= 


cbg 

■+ 



+ 

b( 

s z ") 
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where 


CO = the angular acceleration vector of the spacecraft in 
the coordinates. 


Let the principal moments of inertia of the magnet of the damper be along 
the axes of the coordinate system. Let ^ and 6 ^ be the 

pitch, the roll and the yaw angular displacements of the system 

with respect to the coordinate axes. The relative orientations of 

the /'y'z' coordinates from the coordinates are similar to those 

of the coordinates from the coordinates. Hence, the ly* ^ 

coordinates can he transformed to the x'y'z' coordinates in a manner 
analogous to that given by Eq. (1). Therefore, 


! 




( 8 ) 


Let £0^ be the angular velocity vector of the magnet in the Z 

coordinates. Then and are related by the transformation 






(9) 


For small angular displacements, Eq. (9) can be approximated by 






( 10 ) 
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Hence, Eq. (8) becomes 


3K ^ 




Uj + f T „ “ y 3 ^// ^■*0 


+ ( Jy - I^)(^CO^ U)^ - &,, ^3/^ 






J 



and -Z^ be the principal moments of inertia along the 
coordinate axes. Also, let angular velocity 

vector of the spacecraft expressed in the coordinates. Then 

the equations of motion of the damper magnet are given by 


^et ij, 







+ CT^~-Z^y)^y + 

-h + 


co^- cX>s4>c) - 

b cOy - — 7^ y 

b(a)^ “ ' 


( 12 ) 


where is the angular acceleration vector of the magnet in the 

x'y' s' coordinates and is the Earth's magnetic 

torque vector on the damper magnet. 


It is evident that ^ and ypjj are related by the relation 

-i 

^ ^ ^y)J ^ • 


(13) 


Approximating Eq. (13) by the relation 

= C ec4 , K , ^y)]L 4 , , C 7 + ^a) 


(14) 
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Equation (12) can be rewritten as 



1 





co'g + (r'-jy) Cj;^ 





( 15 ) 


Equations (11) and (15) constitute the approxiniate coupled set of Euler's 
equations for the spacecraft and the datnper magnet. 


LINEARIZED EQUATIONS 

Let a circular orbit be assumed such that 




= the orbital rate = the rate about the ^ -axis 


a a constant 



(16) 


In terms of the Euler angle rates, the angular velocities of the spacecraft 
are then given by 

co^ = 

CO^ = H- Cc^ . 


Let it be assumed that the angles , Oy and 6^^ have a non-zero bias 
such that 

^ 

Oy - -h 


(18) 
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s the bias angles 

= small changes of the attitude angles 
from their bias values. 

and are constants but 

The expressions for the angular velocities as given by Eq. (17) are linearized 
about the bias positions to become 

^rs ~~ 7 

COsj = ^>s ^rs 

(19) 

+ 0r:> ^ys ^ ^ys 

CO J t£ - <^y. S-U. 4- Cir> 6^^ Cas ^ 

~~ "7 CdV S U. 

similarly, the linearized forms of the angular accelerations are given by 

co^ = 

+ ^r.f ( 7 ~ 5^/b CerjBrs " 

COj = Cm ~- 

+ '7 Cm Cm 

+ y,S Im, Cm 7 ^y5 Co3 + 7 ^fs ^ 
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where 

1 ^vs * ^>S 
y . ^:j 

Let it also be assumed that 
is a variable. 



( 20 ) 


^ ^ ‘^y^s 

~ 7 “ 4^s ^ 

+■ 5j C« 6^^ C&3 &Ys “ 7 ^ ^ -f £ U ^Ys, J) 


Similarly, writing 


and 


< = 
< = 


<9, 


^ — 


K. ^ 


( 21 ) 


equations for the damper magnet rates and accelerations analogous to Eqs. 

(19) and (20) are obtained. 

It should be noted that though can be set equal to zero without too 

much error; 0* cannot be. This is because the magnet continuously 
follows the Earth's magnetic field. In a polar or in an equatorial orbit, 
the magnet stays parallel to the field at all times and does not roll. But 
if the orbit Inclination ^ lies between zero and ninety degrees, the mag- 
net roll angle oscillates between O and degrees with twice the orbit- 

al frequency. This is because the Earth's resultant magnetic field makes an 
angle of approximately 2o( degrees with the polar axis at a place with a 
latitude of o( degrees. Thus, can be expressed as 


a 


rs 


~ ^ / - C« 2-y {- ^ 


( 21 ) 


and hence, 


e 


rs 


2^-f) SI 


7 




7 


( 22 ) 
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is the inclination 


where is the time from the ascending node and 

of the orbit. 

Let ^ be the maximum yaw stiffness for an orbiting magnet. Thus, ^ must 
necessarily be the yaw stiffness for a magnet in a polar orbit at the equator. 
Hence, at any arbitrary orbit inclination, , 


T 

= Yaw stiffness = 

7“ 



= Roll stiffness = 


(23) 

TCH 

= Pitch stiffness = 

T“ £ ^ 



where 0( is the local latitude. 

The variation of the absolute value of Cs( over a full orbit can be approx- 
imately written as 



- Crri - 


e 


rs 


(24) 


From these relations, the magnetic torques on the dampers can be expressed 



- j 

it 

- L r c« J 

T 

- L Ts^ ^2 


(25) 


Substituting the expressions for ^ and ^ from Eqns. (19) and (20) to- 
gether with the similar expressions for and Cjp Into Eq. (11), the 
linearized homogeneous form of Eq. (11) can be expressed as 
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( 26 ) 


where 

^/•5 

has been set equal to 

zero. 

The elements 

of the 3x3 matrices 

LriJ 

follows 





s= 

I. 



== 
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'^S,S 2> 

= 



T 

' l,zi 

S£ 
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T 

' 1 , z z 

S.' 

Jy Ccm 


^ 1,23 

a 



^X,3/ 


O 


^J.,32- 

=;i 



~^1, 33 


Coo Co^6>y^ 



= 

t> 




[_ ( ~ Xy J Coo 2. 




^ Cl^ - X,) - 




(27a) 
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2 Z " ~ 7 ^ ^ ^ ^i^s 

^2,23 = (1^- ^ L ^ ^ ^:^s 

^2,3/ “ ~ -^x --^ ) ^ 

^2,22 - ('i’x 

^2^3 ~ f-2'x"-^y}'7,‘^'^ -^-t^ 2 ^ Ctv Sf^ 


(27b) 


4 C^j ile-d SC^&Ys )J 
^3 , /2 ^ ^ 7 )" 

4 3 (Tfrd Ced ^V5 C^'C^ 

— SXw ) + 3 5.. ^ 

C^ ^ys } J 

^ . /3 == ^ 7 / ^rs ^ 2 ^ ^ ^ ^ "" 

- Ce^d ^Ys) ('^ ^ ^/i ^rs 

— ^'ys') 3(^-^ S ^ 
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+ Co-a "*” ^en Xcw. ^ys ^3 

'^3,21 * ^) ( "^^ “ ^^ ) [. 2 SUx. 2-^y^ 

*f 3 ^j>s ^Ys s ^ ^ys X/Ck ^rs)3 

~^2,22 ~ f^*y>) ^ ^ ^ 

— 3 Cea ^ Xjx,^y^ ^S ^S ^ s 

S'Zt^ 

“~ Cn Sa.^ 'S'-tA^ ^Ys) ^ ^Y£ 

+ Ctrj By^ Cfra Sl^ )J 

^3,3/ ^ ('^) '~ ■^y)£i ^ys (27c) 

— Z B^^Ccv ^Ys ^ys ^ ^ ^l>S 

^3, 32 “ ^ ^3 " ^y)£'^‘^ ^ys ^ ^rs ^ ^>5 

“ 3 ^ x>v\^ ^Ks Ccra ^ ^ *5"-tA\. ^ &3 B ^ <S aak^ ^fs 
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T 

V, O 




=r 

o 

zz 

SZ 

- fa 

-Z-3 
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“ •$ •"C'V^ ^erj 

^,3/ 
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o 

32 
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A -T-t4v^ 
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Proceeding similarly, the linearized hcnnogeneous form of Eq. (15) can be 
expressed as 




( 28 ) 


The elements of the 3x3 matrices 
obtained as follows: 


t =- I • -4 , can be 





/2 
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O 

^±. f3 
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’’ S 0y^ 

^ i,Zf 
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^I.IZ 
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ss 

C^6y^ 
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y £ (X^ - Cd-J - Zjt J 



(29a) 
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( -Zy + ~ ^g) ^ ^rs 

s*,3^ = 7 (i^-ij-s^) Si^eis ^ 

^ 2,32 ® 

^2/33 ~ '^C^x~ ^^^ysCn^Ys 

— fijj Sa^ ^ys 

^3.U = (*jtUg-2:y)Ct.^Ks^”^Ks ^ 

%,/2 = f 

+ ^rj -^x y ^irs 

S3,/3 = ^ 

^3,2/ = i ('7/(''^x'-x;;^< 


(29b) 


(29C) 
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= C^f(r:~. 



- c« 

6y^ 

^ 3, 33 

= T S.U, J 








'~^4, 12. 

= 

C- 



~''4, /3 

= 

h> S,.ly^ 



^4,.Z/ 

= 

o 



^4 ,Zx 

— ■ 

— t) CcTJ 



^4,22> 

= 

— i> i -t-w 

Ce-j 


^4,3f 

=• 

o 



^4,3z 


-4xV^' ^ ^ ^ 



^4, 3 3 

— 

— fc> C-erj ^ 

Cej 

^rs 


ifS 


S 2.S 



Csr3 29'^^ 


rs 


(29d) 
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It can easily be verified that assuming a si 05 )le dipole model for the Earth's 
magnetic field, the value of 7“ is given by 


where 


T - 


O. 0/^9 




M 

% 


the strength of the magnet in pole-cms., 
the orbital rates in radlans/sec. 


( 30 ) 


In view of Eqs. (21) and (25), it is evident that the coefficient matrices 

and E^Cj » , of Eqs. (26) and (28) contain time- 

dependent terms. Since these terms are small, and since the time constants 
of the dependent variables are usually of the order of twenty or more orbital 
periods, these time-dependent teims in Eqs. (26) and (28) are replaced by 
their average values over one orbit period. This procedure converts the 
Eqs. (26) and (28) into equations with constant coefficients. 

Let the state vector Ur ^ defined as follows: 

U = ^ (31) 


Then Eqs. (26) and (28) can be combined to become 

EaJ ^ UbJ ^ ^ 


(32) 
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The 


12x12 matrices £AJ and are defined as 


/?J r<^j COJ OOj 

foj LIJ 

[A] = 

CoJ LOJ LT^J Loj 

LoJ [03 LoJ 

£o3 Co] -LiJ 
[oj LoJ Lo] ~[fJ 

(34) 

[T 3 ] [oJ [r^J [rj 
[oJ LS 3 ] LS 4 ] 


where 

QoJ = the null matrix, and 

L X J ® the identity matrix. 
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Let UCj l>e the (12x12) matrix defined by 


Zc3 = - LA . 


( 35 ) 


The time constants and the Garber stability boundaries of the spacecraft- 
magnet system are obtained by confuting the eigenvalues of the matrix 


A-23 



2X2 CYCLES 



















RESULTS 


The results of the linear eigenvalue analysis have been plotted in Figs, A2, 
A3 and A4 and in Figures 4.11, 4.12 and 4.13 which are included 

in the main body of this report. These figures show the effects of the dawp- 
ing coefficient, i , and those of the extra mass added to the satellite 
in its pitch plane to create a difference In the pitch and roll moments of 
inertia of the satellite. 

Figures 4.11, 4,12 and 4.13 show the variation of the pitch damping 

time constant for three different magnet strengths. For a given magnet 
strength ^ M , a given set of the moments of inertia of the satellite and 
for low damping levels, the pitch time constant, as expected, shows a 
linear variation with the damping coefficient when plotted on a log-log 
scale. For large values of the damping coefficient, the magnet Is pulled 
away from the Earth’s magnetic field and the time constant rises. 

Figures A2, A3 and A4 show the corresponding effects of the damping coeffic- 
ient, added mass and the magnet strength on the roll time constant. Though 
the roll time constants are less than the pitch time constants, the effects 
of the roll-yaw coupling between the spacecraft and the magnet is predomin- 
ant here. The linear analysis with zero bias positions shows alternating 
points at vhich roll and yaw instability occurs. But at these unstable 
points, the spacecraft rolls or yaws to a new non-zero bias position which 
Is stable. Thus, the spacecraft will show a small but constant deviation 
from the zero bias positions when the spacecraft is operated at these 
unstable points. 
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SECTION 1.0 


INTRODUCTION 

This report, entitled "Passive Stabilization of the Long Duration Exposure 
Facility" has been, prepared for Langley Research Center under contract num- 
ber NAS 1-13440. It presents the results of a seventeen week study on the 
application of the Magnetically Anchored Rate Damper to gravity gradient 
stabilization of the Long Duration Exposure Facility (LDEF) . The study ob- 
jectives were to perform the analyses and simulations required to investigate 
the use of an existing viscous magnetic rate damper for rate stabilizing the 
LDEF spacecraft. The study itself was broken into three main tasks; Linear 
Performance Estimates, Capture and Damper Requirements, and Performance Pre- 
diction. Each of these tasks was performed for two gravity gradient stabil- 
ization configurations; an axisymmetric configuration for two-axis (pitch and 
roll) stability; and a non-axlsymmetric configuration for three-axis stability. 

The report presents the results by stabilization configuration. 



SECTION 2.0 


SUMMARY AND CONCLUSIONS 

2.1 SUMMARY 

The nominal LDEF configuration and the anticipated orbit parameters are 
shown in Table 2.1. Using these parameters, two linear steady state analy- 
ses were performed; one for two axis (pitch and roll) stabilization, and 
one for three axis (pitch, roll and yaw) stabilization. In each of these 
analyses, the effects of orbit eccentricity, solar pressure, aerodynamic 
pressure, magnetic dipole, and the magnetically anchored rate damper were 
evaluated to determine the configuration sensitivity to variations in these 
parameters. The worst case conditions for steady state were identified, 
and the performance capability calculated. 

Garber instability bounds (a linear instability associated with gravity grad- 
ient stabilized spacecraft) were evaluated for the range of configurations 
and damping coefficients under consideration. 

The transient damping capabilities of the damper were evaluated for both the 
two and three axis configurations, and the time constant as a function of 
damping coefficient, magnet strength, and spacecraft moment of inertia 
determined. The capture capabilities of the danger were calculated, and the 
results combined with the steady state, transient, and Garber Instability 
analyses to select a danqjer design. 

After completion of the linearized analyses, the selected configurations and 
damper design were simulated on a large three axis digital computer program. 
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